(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     14418,        448]
NotebookOptionsPosition[     13761,        429]
NotebookOutlinePosition[     14135,        445]
CellTagsIndexPosition[     14092,        442]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"k", "=", " ", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"g", "^", "2"}], "+", 
     RowBox[{"n", "^", "2"}]}], "]"}]}], ";", 
  RowBox[{"h2", "=", " ", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"k", "[", 
       RowBox[{"g", ",", "n"}], "]"}], "^", "2"}], "+", 
     RowBox[{"1", "i", "*", "w", "*", "u0", "*", "u2", "*", "e2"}]}], "]"}]}],
   ";", 
  RowBox[{"h5", "=", " ", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"k", "[", 
       RowBox[{"g", ",", "n"}], "]"}], "^", "2"}], "+", 
     RowBox[{"1", "i", "*", "w", "*", "u0", "*", "u5", "*", "e5"}]}], "]"}]}],
   ";", 
  RowBox[{"m2", "=", " ", 
   RowBox[{
    RowBox[{"h2", "[", 
     RowBox[{"g", ",", "n"}], "]"}], "/", 
    RowBox[{"[", 
     RowBox[{"u2", "*", 
      RowBox[{"k", "[", 
       RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}]}], ";", 
  RowBox[{"m5", "=", " ", 
   RowBox[{
    RowBox[{"h5", "[", 
     RowBox[{"g", ",", "n"}], "]"}], "/", 
    RowBox[{"[", 
     RowBox[{"u5", "*", 
      RowBox[{"k", "[", 
       RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}]}], ";", 
  RowBox[{"P", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"[", 
      RowBox[{
       RowBox[{"m5", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "+", "1"}], "]"}], "/", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{"m5", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "-", "1"}], "]"}]}], "*", 
    RowBox[{"Exp", "[", 
     RowBox[{"2", "*", 
      RowBox[{"h5", "[", 
       RowBox[{"g", ",", "n"}], "]"}], "*", 
      RowBox[{"[", 
       RowBox[{"d", "+", "t5"}], "]"}]}], "]"}]}]}], ";", 
  RowBox[{"v", "=", " ", 
   RowBox[{"05", "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "+", 
         RowBox[{"m5", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", 
       RowBox[{"Exp", "[", 
        RowBox[{
         RowBox[{"[", 
          RowBox[{
           RowBox[{"h5", "[", 
            RowBox[{"g", ",", "n"}], "]"}], "-", 
           RowBox[{"k", "[", 
            RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", "d"}], "]"}]}], "+", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "-", 
         RowBox[{"m5", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", 
       RowBox[{"P", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "*", 
       RowBox[{"Exp", "[", 
        RowBox[{
         RowBox[{"-", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"h5", "[", 
             RowBox[{"g", ",", "n"}], "]"}], "+", 
            RowBox[{"k", "[", 
             RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}], "*", "d"}], "]"}]}]}],
      "]"}]}]}], ";", 
  RowBox[{"w1", "=", " ", 
   RowBox[{"05", "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "-", 
         RowBox[{"m5", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", 
       RowBox[{"Exp", "[", 
        RowBox[{
         RowBox[{"[", 
          RowBox[{
           RowBox[{"h5", "[", 
            RowBox[{"g", ",", "n"}], "]"}], "+", 
           RowBox[{"k", "[", 
            RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", "d"}], "]"}]}], "+", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "-", 
         RowBox[{"m5", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", 
       RowBox[{"P", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "*", 
       RowBox[{"Exp", "[", 
        RowBox[{
         RowBox[{"[", 
          RowBox[{
           RowBox[{"k", "[", 
            RowBox[{"g", ",", "n"}], "]"}], "-", 
           RowBox[{"h5", "[", 
            RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", "d"}], "]"}]}]}], 
     "]"}]}]}], ";", 
  RowBox[{"x", "=", " ", 
   RowBox[{"05", "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{"[", 
       RowBox[{
        RowBox[{"v", "[", 
         RowBox[{"g", ",", "n"}], "]"}], "+", 
        RowBox[{"w1", "[", 
         RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "+", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "/", 
         RowBox[{"m2", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", 
       RowBox[{"[", 
        RowBox[{
         RowBox[{"v", "[", 
          RowBox[{"g", ",", "n"}], "]"}], "-", 
         RowBox[{"w1", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}]}], "]"}]}]}], ";", 
  RowBox[{"y", "=", " ", 
   RowBox[{"05", "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{"[", 
       RowBox[{
        RowBox[{"v", "[", 
         RowBox[{"g", ",", "n"}], "]"}], "+", 
        RowBox[{"w1", "[", 
         RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "-", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "/", 
         RowBox[{"m2", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", 
       RowBox[{"[", 
        RowBox[{
         RowBox[{"v", "[", 
          RowBox[{"g", ",", "n"}], "]"}], "-", 
         RowBox[{"w1", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}]}], "]"}]}]}], ";", 
  RowBox[{"Ciz", "=", " ", 
   RowBox[{
    RowBox[{"[", 
     RowBox[{
      RowBox[{"-", "2"}], "*", "u0", "*", 
      RowBox[{"k", "[", 
       RowBox[{"g", ",", "n"}], "]"}], "*", 
      RowBox[{"Sin", "[", 
       RowBox[{"a2", "*", "g"}], "]"}], "*", 
      RowBox[{"Sin", "[", 
       RowBox[{"a1", "*", "n"}], "]"}]}], "]"}], "/", 
    RowBox[{"[", 
     RowBox[{"g", "*", "n", "*", 
      RowBox[{"Exp", "[", 
       RowBox[{"S1", "*", 
        RowBox[{"k", "[", 
         RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}], "]"}]}]}], ";", 
  RowBox[{"Cizd", "=", " ", 
   RowBox[{
    RowBox[{"[", 
     RowBox[{
      RowBox[{"-", "2"}], "*", "u0", "*", 
      RowBox[{"k", "[", 
       RowBox[{"g", ",", "n"}], "]"}], "*", 
      RowBox[{"Sin", "[", 
       RowBox[{"a2", "*", "g"}], "]"}], "*", 
      RowBox[{"Sin", "[", 
       RowBox[{"a1", "*", "n"}], "]"}], "*", 
      RowBox[{"Exp", "[", 
       RowBox[{"S1", "*", 
        RowBox[{"k", "[", 
         RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}], "]"}], "/", 
    RowBox[{"[", 
     RowBox[{"g", "*", "n"}], "]"}]}]}], ";", 
  RowBox[{"A", "=", " ", 
   RowBox[{"05", "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "+", 
         RowBox[{"[", 
          RowBox[{"1", "/", 
           RowBox[{"m2", "[", 
            RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}], "]"}], "*", 
       RowBox[{"Ciz", "[", 
        RowBox[{"g", ",", "n"}], "]"}]}], "+", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{
         RowBox[{"[", 
          RowBox[{"1", "/", 
           RowBox[{"m2", "[", 
            RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "-", "1"}], "]"}], "*", 
       RowBox[{"Cizd", "[", 
        RowBox[{"g", ",", "n"}], "]"}]}]}], "]"}]}]}], ";", 
  RowBox[{"B", "=", " ", 
   RowBox[{"05", "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{
       RowBox[{"[", 
        RowBox[{"1", "-", 
         RowBox[{"[", 
          RowBox[{"1", "/", 
           RowBox[{"m2", "[", 
            RowBox[{"g", ",", "n"}], "]"}]}], "]"}]}], "]"}], "*", 
       RowBox[{"Ciz", "[", 
        RowBox[{"g", ",", "n"}], "]"}]}], "-", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{
         RowBox[{"[", 
          RowBox[{"1", "/", 
           RowBox[{"m2", "[", 
            RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "+", "1"}], "]"}], "*", 
       RowBox[{"Cizd", "[", 
        RowBox[{"g", ",", "n"}], "]"}]}]}], "]"}]}]}], ";", 
  RowBox[{"Dd", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"[", 
      RowBox[{
       RowBox[{"m2", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "-", "1"}], "]"}], "/", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{"m2", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "+", "1"}], "]"}]}], "*", 
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", "2"}], "*", 
      RowBox[{"h2", "[", 
       RowBox[{"g", ",", "n"}], "]"}], "*", "t2"}], "]"}]}]}], ";", 
  RowBox[{"Crz", "=", " ", 
   RowBox[{
    RowBox[{"v", "[", 
     RowBox[{"g", ",", "n"}], "]"}], "*", 
    RowBox[{
     RowBox[{"[", 
      RowBox[{
       RowBox[{
        RowBox[{"Dd", "[", 
         RowBox[{"g", ",", "n"}], "]"}], "*", 
        RowBox[{"A", "[", 
         RowBox[{"g", ",", "n"}], "]"}]}], "-", 
       RowBox[{"B", "[", 
        RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "/", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{"y", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "-", 
       RowBox[{
        RowBox[{"Dd", "[", 
         RowBox[{"g", ",", "n"}], "]"}], "*", 
        RowBox[{"x", "[", 
         RowBox[{"g", ",", "n"}], "]"}]}]}], "]"}]}]}]}], ";", 
  RowBox[{"Crzd", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"w1", "[", 
      RowBox[{"g", ",", "n"}], "]"}], "*", 
     RowBox[{
      RowBox[{"[", 
       RowBox[{
        RowBox[{
         RowBox[{"Dd", "[", 
          RowBox[{"g", ",", "n"}], "]"}], "*", 
         RowBox[{"A", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "-", 
        RowBox[{"B", "[", 
         RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "/", 
      RowBox[{"[", 
       RowBox[{
        RowBox[{"y", "[", 
         RowBox[{"g", ",", "n"}], "]"}], "-", 
        RowBox[{
         RowBox[{"Dd", "[", 
          RowBox[{"g", ",", "n"}], "]"}], "*", 
         RowBox[{"x", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}]}], "]"}]}]}], "-", 
    RowBox[{"Cizd", "[", 
     RowBox[{"g", ",", "n"}], "]"}]}]}], ";", 
  RowBox[{"bz", "=", " ", 
   RowBox[{
    RowBox[{"[", 
     RowBox[{"1", "/", 
      RowBox[{"[", 
       RowBox[{"4", "*", 
        RowBox[{"pi", "^", "2"}]}], "]"}]}], "]"}], "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{
       RowBox[{"Crz", "[", 
        RowBox[{"g", ",", "n"}], "]"}], "*", 
       RowBox[{"Exp", "[", 
        RowBox[{
         RowBox[{"k", "[", 
          RowBox[{"g", ",", "n"}], "]"}], "*", "z"}], "]"}]}], "+", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{
         RowBox[{"Cizd", "[", 
          RowBox[{"g", ",", "n"}], "]"}], "+", 
         RowBox[{"Crzd", "[", 
          RowBox[{"g", ",", "n"}], "]"}]}], "]"}], "*", 
       RowBox[{"Exp", "[", 
        RowBox[{
         RowBox[{"-", "1"}], "*", 
         RowBox[{"k", "[", 
          RowBox[{"g", ",", "n"}], "]"}], "*", "z"}], "]"}]}]}], "]"}], "*", 
    RowBox[{"[", 
     RowBox[{
      RowBox[{
       RowBox[{"[", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{"1", "i", "*", 
           RowBox[{"[", 
            RowBox[{"b2", "+", "b2d"}], "]"}], "*", "g"}], "]"}], "-", 
         RowBox[{"Exp", "[", 
          RowBox[{"1", "i", "*", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"-", "b2"}], "+", "b2d"}], "]"}], "*", "g"}], "]"}]}], 
        "]"}], "/", 
       RowBox[{"[", 
        RowBox[{"1", "i", "*", "g"}], "]"}]}], "*", 
      RowBox[{
       RowBox[{"[", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{"1", "i", "*", 
           RowBox[{"[", 
            RowBox[{"b1", "+", "b1d"}], "]"}], "*", "n"}], "]"}], "-", 
         RowBox[{"Exp", "[", 
          RowBox[{"1", "i", "*", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"-", "b1"}], "+", "b1d"}], "]"}], "*", "n"}], "]"}]}], 
        "]"}], "/", 
       RowBox[{"[", 
        RowBox[{"1", "i", "*", "n"}], "]"}]}]}], "]"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.783343686445168*^9, 3.783343742342104*^9}, {
  3.783343814701239*^9, 3.783343828890316*^9}, {3.7833439062613935`*^9, 
  3.783343908019693*^9}, {3.7833446748762007`*^9, 3.783344686540063*^9}},
 FontWeight->"Plain",ExpressionUUID->"501faf64-4892-45fe-8902-f504bc2d8603"],

Cell["\<\
\:8bf4\:660e\:ff1a\:6211\:662f\:901a\:8fc7\:5b57\:7b26\:66ff\:6362\:7684\:529e\
\:6cd5\:5c06\:7a0b\:5e8f\:4e2d\:7edd\:5927\:90e8\:5206\:6539\:597d\:4e86\:ff0c\
\:7136\:800c\:8fd8\:6709\:4e00\:4e9b\:5c0f\:5730\:65b9\:9700\:8981\:4fee\:6539\
\:3002
1. \:62ec\:53f7\:ff1a\:4e2d\:62ec\:53f7\:5728MMA\:4e2d\:4ee3\:8868\:7684\:662f\
\:51fd\:6570\:ff0c\:6bd4\:5982Sin[ x]\:3002\:8be5\:7528\:5c0f\:62ec\:53f7\
\:7684\:5730\:65b9\:662f\:4e0d\:80fd\:7528\:4e2d\:62ec\:53f7\:7684\:ff0c\:6bd4\
\:59821+(x-z)\:4e0d\:80fd\:5199\:6210
1+[x-z]\:3002\:62ec\:53f7\:66ff\:6362\:7684\:5927\:90e8\:5206\:5de5\:4f5c\
\:6211\:5df2\:7ecf\:505a\:597d\:4e86\:ff0c\:5269\:4e0b\:7684\:62ec\:53f7\:8c03\
\:6574\:6211\:76f8\:4fe1\:4f60\:4e5f\:80fd\:505a\:4e86\:3002
2.\:51fd\:6570D\:ff0c\:5728MMA\:4e2d\:5b57\:6bcdD\:5df2\:7ecf\:6709\:5185\
\:7f6e\:5b9a\:4e49\:4e86\:ff0c\:6240\:4ee5\:4f60\:539f\:6765\:7684\:51fd\:6570\
D\:6211\:5168\:90fd\:6539\:540d\:6210Dd\:4e86\:ff0c\:8fd9\:4e2a\:4e0d\:9700\
\:8981\:4f60\:518d\:4fee\:6539\:ff0c\:53ea\:662f\:544a\:8bc9\:4f60\:4e00\:58f0\
\:3002\
\>", "Subsection",
 CellChangeTimes->{{3.7833440802083874`*^9, 3.7833444916358943`*^9}, {
  3.7833447027801647`*^9, 
  3.783344802443229*^9}},ExpressionUUID->"6c67185f-fccb-4b0e-9665-\
2787c42a104c"],

Cell[CellGroupData[{

Cell["", "Subsection",
 CellChangeTimes->{{3.7833444959243956`*^9, 
  3.7833445931426487`*^9}},ExpressionUUID->"2a0ff068-0121-42a5-a132-\
587d7893e7fd"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.7833438857984486`*^9, 
  3.7833439036075277`*^9}},ExpressionUUID->"3f930ffb-9e3f-4fc8-a6b0-\
a99aef4d563b"]
}, Open  ]]
},
WindowSize->{1280, 637},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification:>1.1 Inherited,
FrontEndVersion->"11.3 for Microsoft Windows (64-bit) (March 28, 2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 11583, 371, 451, "Input",ExpressionUUID->"501faf64-4892-45fe-8902-f504bc2d8603"],
Cell[12144, 393, 1265, 21, 176, "Subsection",ExpressionUUID->"6c67185f-fccb-4b0e-9665-2787c42a104c"],
Cell[CellGroupData[{
Cell[13434, 418, 152, 3, 41, "Subsection",ExpressionUUID->"2a0ff068-0121-42a5-a132-587d7893e7fd"],
Cell[13589, 423, 156, 3, 32, "Input",ExpressionUUID->"3f930ffb-9e3f-4fc8-a6b0-a99aef4d563b"]
}, Open  ]]
}
]
*)

